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NOMENCLATURE 


Alphanumeric: 

c. i th chemical species 

H total enthalpy 

N total number of chemical species 

P pressure 

Pr Prandtl number 

r distance from the nozzle center line 

s length along the boundary layer path 

u velocity along s 

v velocity normal to s 

x distance parallel to the nozzle center line 

y distance measured normal to s 

Greek: 

a mixing fraction; =0 for primary flow, =1 for secondary flow 
8 boundary layer thickness 

e eddy viscosity 

i u dynamic viscosity 

p gas density 

a> chemical species production/dissipation term 
Subscripts: 

b£ boundary layer 

core primary flow 

e boundary layer edge 

gg gas generator exhaust (secondary flow) 

inj injectant (secondary flow) 

trans transpired fluid (secondary flow) 

T turbulent 

w wall 

Superscripts: 

j j=0 for planar symmetry and j-1 for axisymmetric 
fluctuating with time 
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NOMENCLATURE (concluded) 


Abbreviations: 

DROPMEX 

MABL 

MABL-E 

MABL-F 

MABL-K 

P.D.E. 

PNS 

SCAP 

SOW 

SSME 

STME 

TDE 

TDK 

VIPER 


computer code for injector performance prediction 
Mass Addition Boundary Layer module of TDK 
MABL with equilibrium chemistry 

MABL with chemical compositions frozen at the stagnation value 

MABL with finite rate chemical kinetics 

partial differential equations 

Parabolized Navier Stokes 

Spray Combustion Analysis Program 

Statement of Work 

Space Shuttle Main Engine 

Space Transportation Main Engine 

TDK with equilibrium chemistry 

Two Dimensional nozzle analysis code with chemical kinetics 
Viscous Interaction Performance Evaluation Routine computer code 
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1.0 INTRODUCTION 


An advanced version of the Two-Dimensional Kinetics (TDK) 1 computer program was 
developed under contract NAS8-36863 and released to the propulsion community in early 
1989. Exposure of the code to this community indicated a need for improvements in certain 
areas. 


In particular, the TDK code needed to be adapted to the special requirements imposed 
by the Space Transportation Main Engine (STME) development program. This engine utilizes 
injection of the gas generator exhaust into the primary nozzle by means of a set of slots. The 
subsequent mixing of this secondary stream with the primary stream with finite rate chemical 
reaction can have a major impact on the engine performance and the thermal protection of the 
nozzle wall. In attempting to calculate this reacting boundary layer problem, the Mass 
Addition Boundary Layer (MABL) module of TDK was found to be deficient in several 
respects. For example, when finite rate chemistry was used to determine gas properties, 
(MABL-K option) the program run times became excessive because extremely small step sizes 
were required to maintain numerical stability. A robust solution algorithm was required so 
that the MABL-K option could be viable as a rocket propulsion industry design tool. Solving 
this problem was a primary goal of the Phase I work effort under Contract NAS8-39048. 
Phase I was conducted between 1 August 1991 and 31 December 1991. The Statement of 
Work (SOW) for Phase I is shown in Table 1. 

Towards the end of the Phase I effort it became apparent that much labor could be 
saved if the TDK code was made more compatible with the analysis tools being used for 
calculating the STME power cycle, and for calculating distributed energy release efficiency. 
For this reason a Phase II effort was initiated when funds became available in September 1992. 
Phase II of the contract was conducted between 16 September 1992 and 19 March 1993. The 
SOW for Phase II is shown in Table 2. 


1. Nickerson, G.R., Coats, D.E., Dang, A.L., Dunn, S.S., and Kehtarnavaz, H„ Two- 
Dimensional Kinetics (TDK) Nozzle Performance Computer Program, Volume I, II, and 
III, Software and Engineering Associates, Inc., NAS8-36863, dated 31 March 1989, 
prepared for George C. Marshall Space Flight Center. 
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Table 1. Phase I Statement of Work, NAS8-39048 


STATEMENT OF WORK 


The Two-Dimensional Kinetics (TDK) computer program is a large 
software package oriented to predict the performance of a liquid 
rocket thrust chamber. Many options are available for a variety of 
problems, and several boundary layer modules are coupled with the 
code. Substantial modifications to the program have been completed 
recently, and a rigorous documentation set has been prepared to use 
this analytical capability for the development of projected future 
engine designs. 

The contemplated effort is directed to improve some prediction 
features, supporting associated design elements in the thrust 
chamber, to provide support for the code operation on different 
computer systems, and to resolve difficulties occurring during 
program execution. The individual tasks are outlined below: 

Tasks 

1. MABL Restart 

Add a restart capability to the Mass Addition Boundary Layer 
(MABL) module so that after TDK has been run, MABL can be 
executed from a saved file using different options. 

2. MABL-K Stability 

Develop a robust kinetic algorithm for MABL by applying the 
method used to make TDK numerically reliable at high chamber 
pressures. 

3. MABL-K with Injection 

Make MABL with kinetics applicable to nozzles with tangential 
slot injection and transpiration cooling. 

4 . Technical Support and Documentation 

Technical support of the code will consist of the following 
activities: 

o Error correction 

o Hot-line response to user problems 

o Program modifications that may be necessary to make the code 
fully useful 

o Updates to the input preprocessor and postprocessor 
o Additions and improvements to the TDK program documentation, 
version March 1989. 
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Table 2. Phase II Statement of Work, NAS8-39048 


STATEMENT OF WORK 


The Two-Dimensional Kinetics (TDK) computer program is a large 
software package oriented to predict the performance of a liquid 
rocket thrust chamber. Many options are available for a variety of 
problems, and several boundary layer modules are coupled with the 
code. Substantial modifications to the program have been completed 
recently, and a rigorous documentation set has been prepared to use 
this analytical capability for the development of projected future 
engine designs. 

The projected six months effort will support an optimum. STME engine 
design based on the interaction of the performance prediction codes 
with the power cycle and energy release programs. At present, only 
the thrust chamber performance is determined with the TDK/MABL code 
based on provided interface conditions. However, engine component 
optimization depends strongly on the flow distribution in the 
selected gas generator cycle approach. To accomplish this goal, the 
interaction with 'Power Cycle’ and energy release codes must be 
coordinated. Especially, losses from the injection element 
operation and the boundary layer behavior with flow injection are 
of interest. The individual tasks are outlined below: 

Tasks 

1 Provide technical guidance for the program operation (TDK, 
MABL, BLM, Preprocessor, Processing of Results) . 

2 Calculate the mass flow in the boundary layer and the 
associated profiles, including injected mass flows. 

3 Activate the shock calculation with the tri-propellant option. 

4 Prepare information from TDK program solutions for use in an 
engine power cycle code (Efficiencies for C and ISP, Mass 
Fractions, etc. ) 

5 Support the effort to standardize the energy release method 
and to implement the respective features into the TDK program. 
This includes review of the DROPMIX program principles. 

6 Prepare a file holding x-y arrays of significant parameters 
for result visibility, plotting, and optional printout for TDK 
results along the wall, the center line, the boundary layer 
edge, profiles normal to the wall, and in the exit plane. 

7 Update the preprocessor program, MSFC version, to include all 
input items related to the latest TDK program changes. 

8 Update the existing TDK documentation as required. 

9 Quarterly progress reporting and a final report is required. 
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2.0 SUMMARY OF RESULTS, PHASE I 


MABL-K Analysis 

The boundary layer equations for compressible turbulent flow including the effects of 
finite rate chemical kinetics can be derived from the time dependent Navier-Stokes equations 
using the Reynolds time-averaging procedure and the usual boundary layer order of magnitude 
assumptions. As presented below, the boundary layer equations are written in a curvilinear 
coordinate system in which s is the wetted length along the wall and y is measured normal to 
it (x is axial distance measured along the centerline). A bar over a quantity denotes that it is a 
time averaged quantity. It has been assumed that the lateral and transverse curvature terms 
can be neglected. The conservation equations are: 


Continuity 


I (pK ) + i[ ( p v + p' w 'K 


= 0 


( 1 ) 


Momentum 


P u ^+ (pv +pV) 




dn 


( 2 ) 


Energ y 


„ dH ^ , n , -TTs <?H d 
pu^ + cpv + p v 


Pr Pr^ 


*1 

37 



' 1 


' , 1 

■ 

du 

p 

1 " Pr 

+ pe 

1 " 

k T J 

. 

u ^. 


(3) 


Species Conservation 


dc 


dc [ _ d 


pu^ + (pv + pV)^=^ 


Pr + p£ Pr 


k]$b* 


i = 1,2,...N 


(4) 
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where j=0 for 2-D planar flow, j=l for axi symmetric flow. The quantities ju and Pr, are the 
dynamic viscosity, and Prandtl number, respectively. The quantities e and Pr T are their 

turbulent counterparts. 

Boundary conditions for these P.D.E.'s are applied at the boundary layer edge, and at the 
wall, as described below. 

At the edge , y = <5 

Flow conditions here are inviscid and the boundary conditions are obtained using the 
TDK code, i.e. 

u = u 

e 

H = H 

e 

c. = c. 

1 1 

e 

In practice the pressure distribution, rather than the edge velocity, is specified and the velocity 
distribution is found from the pressure gradient by integrating Bernoulli's equation. 

At the wall , y = 0, the boundary conditions are: 
u =0 

-r- = 0 , for adiabatic wall, or 

o y 

T = T w , for specified wall temperature distribution. 


The effects of transpiration cooling are approximated by modifying the wall boundary 
condition. At the wall 

u I = 0 , el = 0 , and ft). I = small 
1 w 1 w i 1 w 

so that equation (4) becomes 


(wall side) 
(9c i 


(p v) 


W 


W 


Thus 


(p v) dc. = 
r w 1 


(b.l. 

d 

Pr <9y 




side) 

dc j 

L w. 

<9q 

w 


T930345 


5 



The species concentrations transpired from the wall are c. . At the boundary layer side of 
the wall the species concentrations are c. . Integration across the wall boundary element 
gives 


(c ; 


i,w 


c. , ) (p v) 

i,trans w 


_ [ J_ dc \ 
Pr W 


so that 


<9c i 

w 



' Pr' 

w 



(p v) (c. - c. ) 

r w i,w l.trans 


' w 


(5) 


Equation (5) is the boundary condition that is applied for a transpiration cooled wall. The 
term (p v) is the rate at which mass is transpired across a unit area at the boundary, and this 

value must not be too large if the method is to work properly. It must be input for each 
species at each station along the transpiration cooled sections of the wall, i.e. 

N 

(p v) = X (p v) c. 
r w . . r 'w l.trans 
i=l 


Nozzles that utilize a tangential slot injection scheme can be treated by the MABL 
analysis if the following assumptions are made: 

1 . The injectant must be introduced parallel to the nozzle wall (tangential 
injection). The velocity of the injectant cannot exceed the velocity of the core 
flow at the point of injection. 

2. The pressure of the injected gas must match that of the core flow. 

3. The mass flow rate of the injected gas is small (less than ~5%) compared to the 
total engine mass flow rate. 

These assumptions are required for the analysis if it is to be consistent with the thin 
shear layer approximations. The application of the method is fully described in Appendix B of 
Reference 1, Volume III. 

The difficulty in integrating the kinetic boundary layer equations (1) through (4) results 
from the source term o>, appearing in equation (4). This term can be extremely sensitive to the 
chemical species concentrations, q, leading to numerically "stiff" behavior. Thus, an 

"implicit" PDE integration scheme with some type of active step size control is required. 
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The method used by MABL-K is first order simple back difference. The method 
requires the inversion of a large matrix (see Reference 1, Vol. I, Section 6). The coefficients 
in this matrix include terms of the form 


d C0j / d Ci 


and these are found analytically, rather than by numerical differencing. Stabilizing the 
MABL-K algorithm required that the implicit method be rederived to assure its correctness, 
since any error in the formulation could adversely affect numerical stability. Furthermore, 
cases were run where the derivative expressions were evaluated by numerical differencing. 
The matrix coefficients obtained in this manner were compared to the matrix coefficients 
actually used by the code, i.e. those obtained by formal differentiation. 


The MABL-K was reworked from end-to-end in an effort to improve its numerical 
stability. However, the most fruitful result was achieved by bounding the net change in 
species concentrations during each integration step. Species integration is of the form 


where 


c. i 
i,n+l 


c. + Ac^ 


,< 1 ) 


i ,n+l 


I C. =1 
i 

1A c. = 0 . 

i 


The predicted Acj^ . were bounded so that no individual change exceeded a few 

percent per step. If one species is affected, then all are affected, since atomic conservation 
must be maintained. 


The modifications discussed above have decreased run times required for the MABL-K 
module by factors ranging from 2:1 to 10:1, while increasing both the accuracy and the relia- 
bility of the calculations. 

Figures 1 through 4 show example results obtained using the MABL-K/TDK code to 
calculate boundary layer properties for the SSME (see Reference 1, Section 11.4 for a 
discussion of performance calculations for this engine). For the purposes of these calculations 
an adiabatic wall has been assumed. Edge conditions for the boundary layer were provided by 
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a two dimensional kinetics (TDK) inviscid calculation. MABL-F, MABL-K, and MABL-E 
calculations were then carried out and the results compared (the restart option was used; see 
Table 1, Task 1). The assumptions that differ between the calculations are: 


MABL-F gas in the boundary layer is frozen at the stagnation composition 

MABL-K gas composition is determined by finite rate chemistry (consistent with 

TDK) 

MABL-E the boundary layer is everywhere in a state of chemical equilibrium (gas 
composition is not required). 

In Figure 1 it can be seen that the adiabatic wall temperature as determined by 
MABL-K does not track the equilibrium (MABL-E) values. A similar result is found for 
displacement thickness (see Figure 2). However, the MABL-K predictions for momentum 
thickness and integrated shear stress do track the equilibrium (MABL-E) values (see Figures 3 
and 4). For this case it can be concluded that if wall heat transfer is important, then a 
MABL-K calculation may be required. However, a MABL-E calculation is adequate for 
performance prediction. 


The MABL-K calculation shown above required less than a minute on an IBM 
RISC/6000 320H machine. The calculation would have taken over eight minutes on the 
previous version of the code. 

The MABL-K result for adiabatic wall temperature shown in Figure 1 has been verified 
by repeating this calculation using the VIPER code 2 , which provides a PNS computation for 
the flow field downstream of the nozzle throat. 


2. Kawasaki, A.H., Berker, D.R., Coats, D.E., Dunn, S.S., and Nickerson, G.R., Viscous 
Interaction Performance Evaluation Routine for Two-Phase Nozzle Flows With Finite 
Rate Chemistry, VIPER 2.0, Computer User’s Manual, Software & Engineering Asso- 
ciates, Inc., PL-TR-92-3053, January 1993 prepared for Phillips Laboratory. 
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Figure 1. Adiabatic Wall Temperature vs Length. SSME 
(Note: all calculations assume the same edge velocity.) 


MABL-r 

MABL-K 



X (IN) 


Figure 2. Displacement Thickness, 8 *, vs Length. SSME 
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IOOO. 7000 - 3000. 4000. 61)1)0. 6000. , -Unn .020 .040 .060 -080 .100 


MABL-E 

MABL-K 



Figure 3. Momentum Thickness, 0, vs Length. SSME 



°-20 ■ 0 .0 20.0 40.0 60.0 80.0 100-0 3 20 . G 

X (IN) 


Figure 4. Integrated Wall Shear Stress, t, vs Length. SSME 
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3.0 SUMMARY OF RESULTS, PHASE II 


A new Users Manual 3 has been prepared that includes modifications to the TDK/MABL 
code that were carried out under this contract. This document replaces Reference 1, Volume 
III. 


In an effort to coordinate with other computer codes being used to evaluate power cycle 
behavior, and energy release effects, several auxiliary calculations were added to the 
TDK/MABL code. Tasks 2 and 4 were included in the SOW for this purpose. The methods 
used are described below. 

In order to calculate the mass flow in the boundary layer and the associated profiles, 
including injected mass flows, the "mass flow in the boundary" layer was taken as: 

5 

m^= j pu2 7rr e dy (6) 

o 

where 8 is that point in the boundary layer where the velocity attains a value that is 99.5% of 
the edge velocity, U , i.e., 

y = 5 when (u/U ) = .995 (7) 

The above definition is arbitrary since the boundary layer profiles attain the edge values 
asymptotically. 


3. Nickerson, G.R., Berker, D.R., Coats, D.E., and Dunn, S.S., Two-Dimensional Kinetics 
(TDK) Nozzle Performance Computer Program Volume III, Users Manual, Software 
and Engineering Associates, Inc., NAS8-39048, dated 31 March 1993, prepared for 
George C. Marshall Space Flight Center. 
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If the nozzle is designed with injection ports along the wall, then for binary mixtures 
(MABL-F and MABL-E) equation (4) is replaced by the following equation: 

Mixture Conservation 


„ da , , , — Tn da d 

puj-tfpv+p v 



( 8 ) 


Equation (6) can then be written as: 


5 

rh^ = J a p u 2 ^ r £ dy + 


5 

\ (1-a) p u 2 K r dy 
o ° 


(9) 


Equation (9), above, separates the boundary layer flow into two components; injectant 
(o^l) and core flow (a=0). It is also required that a—0 at y —8 because the boundary layer 
must contain all of the injectant. 

Thus, the main flow in the boundary layer can contain two components: 


^gg.b^ + ^core.t 


When the injection occurs at a discrete position along the nozzle wall, X. , as is the 

mj 

case with tangential slot injection, then 


= 0 


X<X. 

mj 


m 


gg>b^ 


= m. 


mj 


that is 


8 

jap u 2 7T r dy = m. . 
r e J inj 


X>X. 

mj 


X <X. 
inj 


( 10 ) 


which is a constant value. 
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The following plot files containing boundary layer profiles are provided at up to 20 
specified locations along the nozzle wall (see Reference 3, Vol. II, Section 9.9). 


\l 

m 

gg’bz 


vs y (for MABL-F, MABL-K, and MABL-E options using 

equation (6)) 

vs y (for MABL-F and MABL-E options using equation (10)) 


where 

0 < y < <5. 

Since a pertains only to binary mixtures (MABL-E and MABL-F options), only the m^ 
vs y profiles are provided for kinetic boundary layer (MABL-K) calculations. 

A plot file is also provided that gives 
vs x - 

Application to the STME 

During the course of the work effort, the TDK/MABL code was used to obtain 
performance calculations for a number of rocket engines. Four engines in particular were 
analyzed: 1) the SSME, 2) the P&W 40K thrust level STME test engine, 3) the 583K STME 
design, and 4) the 650K STME design. Analysis results for this latter engine are presented 
here to demonstrate this newest version of the TDK/MABL code. 

The STME design utilizes a combustion chamber extending to a supersonic expansion 
ratio of 7:1 where an expansion nozzle is attached. The configuration is shown in Figures 5a 
and 5b. The nozzle is cooled by gH 2 flowing inside of coolant tubes forming the wall contour. 
Additional cooling is provided by injection of gas generator gases, consisting of primarily gH 2 
from the low mixture ratio reaction. The joint separating the chamber and nozzle is cooled 
with a small amount of gas generator turbine exhaust that is injected normal to the primary 
flow without being choked (subsonic film). The remainer of the turbine exhaust is injected 
parallel to the primary flow as a supersonic film. The mass flow of the supersonic film is 
approximately five times greater than the subsonic film mass flow. The nozzle injector region 
is contained in the dashed circle shown in Figure 5a. An enlarged view of this region is shown 
in Figure 5b. 
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Figure 5a. STME Chamber/Nozzle Configuration 



Convective 
Coolant Passage 

gH 2 


Figure 5b. Detail A: STME Nozzle Injector Region 
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Input data defining the STME geometry and operating conditions is given in Table 3. 
The analysis consists of a two dimensional calculation for the thrust chamber core flow, 
followed by an analysis of the wall boundary layer flow. The TDK code treats the core flow 
as an inviscid chemically reacting process. A wall boundary layer is calculated by the MABL 
code using edge conditions supplied by TDK. Since the flow is found to be in a state very 
near to chemical equilibrium, equilibrium chemistry is assumed; i.e., TDE and MABL-E 
results were obtained and are presented here. 

At the STME chamber/nozzle interface a small amount of turbine exhaust is injected 
normal to the wall, followed approximately three inches downstream by tangential injection of 
a larger amount. Slot injection effects are calculated by the boundary layer analysis (MABL). 
Both slots are treated as tangential injection, since the code has no normal injection capability 
other than injection by transpiration. Thus, flow injected at the first slot is assumed to be 
turned parallel to the core flow. Since the injection flow rate at the first slot is low, and occurs 
behind a rearward facing step, this assumption should give a reasonable approximation of the 
boundary layer properties. 

The TDK/MABL computational procedure requires that the pressure of the injectant be 
equal to the boundary layer edge pressure at the injection position. Thus, the procedure 
calculates the injection nozzle expansion ratio that provides the matched edge pressure. The 
flow from both slots have the same stagnation temperature (see Table 3, T0GG=2* 1 224 °R in 
$ODE), but have different stagnation pressures (P0GG=77.,204. psi for slots 1 and 2, respec- 
tively). The flow is assumed to expand from these values isentropically to match the pressure 
of the adjacent core flow as calculated by the TDE analysis. Computer print-out showing the 
results of these expansion calculations is presented in Table 4. Included in this output is the 
ratio of the mass flow per unit area of the injectant to the mass flow per unit area of the 
boundary layer edge flow. At the first slot the coolant is injected normal to the edge flow. 
However, it can be seen that the ratio of injectant to edge mass flow per unit area is only .166, 
which implies that flow turing effects can be ignored (the analysis requires this). 

Table 4 also shows that the edge pressure at injection is 67.26 and 56.75 psi for the first 
and second injectants, respectively. This give Mach number values after expansion of .447 
and 1.489. Hence, the first slot is said to have "subsonic" injection. The second slot expands 
the film to approximately Mach 1.5, which is believed to provide sufficient margin so that the 
injectant ports will remain choked. Since it turns out that the primary nozzle mass flow is 
predicted to be 1482 lbm/sec, the slot mass flows are 7 and 35 lbm/sec, respectively (see Table 
3, input item FXMGG=.004697, .023687 in $MABL). 
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Table 3* TDE/MABL-E Input Data for the STME Design 
at the 650K Thrust Level 


TITLE F/S NLS 650K ENGINE: GG SUBSONIC AND SUPERSONIC INJECTION 
DATA 
SDATA 


ODE = 1, ODK = 0, TDE = 1 , HABL = 1, IMABL = 3, MABIE 
ASUP = 7. ,43. 09 7, NASUP=2, 

RSI = 6.9, 

X I C< 1 ) = 5.36, NXIC = 1, 

RWTU=0 .494, RWTD= .2, RI = 1.598, 

THETA= 29.228, THE=8.312, THETAI=25 .42 , 

I WALL=4, I TYPE=1 , 


ECRAT =2 .682, A$UB( 1 )=2. 682, NASUB = 1, 

RS(2) = 1.252232, 1.541037, 1.888917, 2.267946, 

2.661167, 3.058127, 3.452169, 3.838668, 
4.214417, 4.576776, 4.924005, 5.254459, 
5.567102, 5.860643, 6.134528, 6.388157, 
6.564863, 

ZS ( 2 > = 0.498084, 0.973261, 1.547740, 2.203318, 

2.935726, 3.742109, 4.620865, 5.571712, 
6.592352, 7.683121, 8.841478,10.067299, 
11.358459,12.715493,14.135952,15.618718, 
16.777452, 


NWS = 18, 


XWGG=2 .836, 3.267, NWGG=2, AFWGG=2*1, EQLGG=T , 
N20NES = 2, XM(1 ) = .925, .075, 

TR I =T , OFCORE= .87827, .83333, 

FFCORE=. 12173, .16667, 

GFCORE=2*0, 

OFGG = .47617, FFGG = 0, GFGG = .52383, 


SEND 

REACTANTS 

O 2. 

H 2. 

H 2. 


100 . 

100 . 

100 . 


= T, 


-2899.1 
-1963. L 
-1963. L 


..2 ZONE 


100.330 
40. 33 F 
40.33G 


NAMELISTS 
$00 E 

RKT = T, 

P = 2250., XP(1) = 2*1, PSIA = T, 

P0GG= 77., 204., T0GG=2* 1 224 . , DELH=2*52 . 265 , 

SEND 

STRANS 

MP = 150, 

DRMIN = .0005, 

SEND 

SMOC 

EXITPL = T, 

NC = -1, 

SEND 

SMABL 


DX I =2E -4 , NDXI=50, NYI=115, DXLIM(I) = .002, .02, 

LDXL I M=2, XLIM=5, 10, 

FXMGG=. 004697, .023687, 

ADBAT C=0, 

OFC=2, D I $TRB=0, 

XC0=-2. 58986, 

XCE=2 . 836, 

ETAC=1 , 

NTQW=25 , 

XTQW= -1000., -2. 6, -1.8, -1.4, -1.1, -0.6, -0.2, -0.1, 

.05, 0.23, 0.85, 1.3, 2., 2. 3, 2. 75, 3. 24, 4. 35, 5. 84, 
7.33,8.82,10.31,11.79,13.28,14.77,15.98, 

TQW = 1100. ,1100. ,1350. ,1300. ,1325. ,1250. ,1350. ,1300., 

950.. 650.. 850.. 800.. 650.. 600.. 1200.. 1215.. 1197.. 1523 

1684.. 1759.. 1793.. 1806.. 1811.. 1803.. 1744., 

ITHKBL = 1, 

SEND 


1.149 

.0709 

.0709 
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Table 4. TDE Print-out of Injection Conditions 


GG MULT I -SLOT OPTION FOR UPPER WALL 


SLOT 

XGG 

(X/R* ) 

RGG 

(R/RM 

POGG 

(PSD 

T0GG 
(DEG R) 

PGG 

(PSD 

TGG 
(DEG R) 

RHOGG 

(LBM/FT3) 

VGG 

(FT/SEC) 

(RHOV)GG 
/ (RHOV)E 

MGG 

AE/AT 

ISP/GG 

(SEC) 

1 SP/EXIT 
(SEC) 

1 

2.836 

2.610 

77.00 

1224. 

67.26 

1179. 

2 .045E • 02 

2049. 

.166 

.447 

1.459 

.00 

264.39 

2 

3.267 

2.829 

204.00 

1224. 

56.75 

859. 

2.368E-02 

5839. 

.617 

1.489 

1.171 

240.59 

281.03 


The wall boundary layer, including interaction of the injectant streams with the primary 
flow, is calculated by the Mass Addition Boundary Layer (MABL) code. Some results 
obtained from these calculations are presented in Figures 6 through 14. The plotted results 
shown in these figures are of two types; 1) boundary layer profiles at a given axial locations, 
and 2) boundary layer parameters such as edge values, wall heat flux, wall shear stress, body 
displacement, etc. vs axial location. 

Figure 6 shows boundary layer profiles for mass flow, a, u/U , and temperature at an 

axial position immediately behind the first slot. The approach boundary layer has been 
displaced, and a slug of injectant has been inserted between it and the wall, as can be seen in 
the plot of y/y^ vs a. The plots of y/y e vs rh^ and rh gg are obtained using equations (6) and 

(10). The vertical line terminating the m ^ profile shows that 7 lbm/sec of turbine exhaust has 

been injected. 

Figure 7 shows these same plots at an axial position immediately behind the second slot. 
The plot of y/y e vs a shows a mixing profile for the first injectant that is superimposed on a 

new slug of injectant from the second slot. The plot of y/y^ vs rn^ shows that a total of 42 

lbm/sec of turbine exhaust has been injected (7 + 35). 

Figure 8 shows these same profiles at the nozzle exit, which is 7.8 feet from slot 
number 2. The plot of y/y vs a shows that at the wall (y/y = 0) approximately 30% of the 

C 6 

mass flow is injectant and 70% is entrained edge flow (of course these are reacted to a state of 

equilibrium since the MABL-E calculation is used here). The calculation indicates that the 

film stays within four inches of the wall (y = 3.9 inches). 

6 
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Figures 9 through 14 are plots of selected boundary layer parameters vs axial position. 
Figure 9 gives m ^ vs X, and can be compared to the profiles given for this quantity in Figures 

6, 7, and 8. The injection positions are seen in Figure 9. Note that at the nozzle exit the 
boundary layer has entrained 240 out of 1482 lbm/sec, or 16% of the total flow. (42 lbm/sec 
of this is injectant). 

Figures 10 through 14 give the user specified wall temperature, integrated wall surface 
heat flux, surface integrated wall shear stress, boundary layer thrust decrement, and body 
displacement vs axial position. The nozzle wall temperatures shown in Figure 10 should be 
regarded as preliminary, and were obtained by an independent analysis provided to SEA, Inc. 
by the NASA/MSFC. This is the profile listed in the TDK input. Table 3, arrays XTQW and 
TQW. Figure 1 1 shows the predicted wall heat pick-up obtained using the wall temperatures 
shown in Figure 10. It should be noted that wall temperature variations of a few hundred °F 
have only a mild effect on this quantity because the adiabatic wall temperature is so high (heat 
flux is roughly proportional to the difference between the adiabatic wall temperature and the 
input wall temperature). It can be seen that injection of the relatively cool turbine exhaust is 
predicted to give a large reduction in wall heating. 

Figure 12 shows the surface integrated wall shear stress. Injection of the relatively low 
velocity turbine exhaust is predicted to give a large reduction in this quantity also, but not so 
large as the reduction in wall heating. A reduction in wall shear stress provides an increase in 
thrust, which in this case could be as much as 5000 lbf. 

Figure 13 shows the predicted boundary layer loss. In this figure the stream thrust of 
the primary injectant can be seen as a vertical line at the injection position (X = 22.5 inches), 
and provides a thrust of - 7500 lbf (this is the force component aligned with the nozzle axis). 
As the injectant expands within the primary flow, it provides additional thrust which fully 
off-sets drag loss until the X ~ 75 inch position is reached. Because of the injectant, the 
boundary layer provides a net gain in engine thrust of about 7000 lbf. 

Figure 14 shows body displacement. It can be seen that one inch of displacement is 
sufficient to accommodate the boundary layer and the injectants. Note that the displacement 
thickness at the nozzle throat is negative due to the large amount of wall cooling taking place 
there, so that the gas density in the boundary layer is increased significantly. The "body 
displacement" is required to adjust the potential to the real wall contour. 
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Figure 9. m . vs X for the 650K STME 




Figure 10. Input Wall Temperature vs X. 650K STME 






Figure 11. Heat Flux Integrated Over the Wall Surface Area vs X. 650K STME 




Figure 12. Wall Shear Stress Integrated Over the Wall Surface Area vs X. 650K STME 






Figure 13. Boundary Layer Thrust Decrement vs X. 650 STME 






Figure 14. Body Displacement vs X. 650K STME 



Figures 15 through 20 show TDK calculations of pressure and temperature along the 
wall (which is the edge condition for the boundary layer) and the flow centerline. In Figure 15 
and 16 a slight compression can be seen at the end of the circular arc used to define the 
downstream side of the throat. This disturbance travels to the flow centerline where it is 
reflected. It can be seen at the X/r*=ll position in Figures 17 and 18. As would be expected 
in an axially symmetric flow, the disturbance is stronger when it reaches the flow center. 

Figures 19 and 20 show pressure and temperature across the nozzle exit plane. The 
horizontal line at the top of Figure 20 represents the interface (or slipline) separating the core 
flow into inner and outer regions. The outer region contains 7.5% of the flow, and this flow is 
contained within the small annulus shown in Figure 20. 

The results plotted in Figure 6 through 20 discussed above were obtained using the x-y 
arrays requested in Task 6 of the Phase II SOW. 

Task 4 of the SOW specified that certain output values be printed for use in an engine 
power cycle code. Example output for the 650K STME calculation is shown in Table 5. Note 
that the primary flow has been separated into two zones 1) an outer region containing 7.5% of 
the flow at a mixture ratio of 7.2149:1, and 2) an inner region containing 92.5% of the flow at 
a mixture ratio of 5:1. The combined over-all mixture ratio for the core flow is 7:1. Table 5 
gives results based on the assumption of one-dimensional isentropic flow. 


Table 5. Performance Summary for the 650K STME 
(one-dimensional, core flow) 






ODE 

ISP & CSTAR SUMMARY 



ZONE 

ODE ISP 

OOF ISP 

ODE 

CSTAR 

ODF CSTAR ZONE MASS FLOW 

M.R. 

OX. FRAC 

1 

2 

449.186 

460.286 

421.249 

449.064 


7353.3 

7870.4 

7217.8 .92500 

7776.1 .07500 

7.2149 

4.9999 

.8783 

.8333 


ODE MASS AVERAGED 
ODF MASS AVERAGED 
ODE MASS AVERAGED 
ODF MASS AVERAGED 

ISP 

ISP 

CSTAR = 
CSTAR = 

450.019 

423.335 

7392.1 

7259.7 






MASS AVERAGED 0/F 

= 6.99357 
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Figure 16. Temperature (°R) vs X/r* Along the Boundary Layer Edge. 650K STME 




Figure 17. Log of Pressure (psi) vs X/r* Along the Centerline. 650K STME 





Figure 18. Temperature (°R) vs X/r* Along the Centerline. 650K STME 






HNOZ 



Figure 20. Temperature (°R) vs r/r* at Nozzle Exit. 650K STME 




The overall performance summary for the 650K STME is given in Table 6. These 
results are intended only as an example. From this information the kinetic, two-dimensional, 
and boundary layer losses and their associated efficiencies with respect to ODE solutions can 
be identified. The energy release loss, consisting of the microscale mixing and evaporation 
processes, are not included in this analysis. A planned JANNAF workshop for the "Standard- 
ization of Injector Performance", which was approved by the Combustion Subcommittee, had 
to be canceled after considerable preparation by several parties involved, since the workshop 
coordinator did not receive the proper support from his organization. The requested review of 
the principles, used in the DROPMIX computer program 4 , for discussion at the workshop is 
presented in Appendix A. One method, concerning the energy release, has been provided by 
means of the SCAP module coupling with the TDK program under an SBIR phase II activity 5 
for hydrocarbon propellants. This information is available to the government only at this time. 


4. Nurick, W.H., DROPMIX-A PC Based Program for Rocket Engine Injector Design, 
27th JANNAF Combustion Meeting, Vol. n, pp 435-468, Cheyenne, Wyoming, 
November 1990. 

5. Nickerson, G.R., and Johnson, C.W., The Chemical Kinetics of LOX/Hydrocarbon 
Combustion, SBIR 88-1, Phase II, Contract NAS8-38468, Software & Engineering 
Associates, Inc., Carson City, Nevada, 25 July 1992. 
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Table 6. Performance Summary, 650K STME 


PRIMARY CHAMBER OPERATING CONDITIONS 


CHAMBER PRESS 

[PSIA] 

2250.000 

CHAMBER TEMP 

[R] 

6659.095 

MIXTURE RATIO 

M 

6.993572 

H (OXID) 

[CAL/MOLE] 

-2899.000 

H (FUEL) 

[CAL/MOLE] 

-1963.000 

HCHAM (ODE) 

[ BTU/LB ] 

-309.6788 

DELH (AVERAGE) 

[ BTU/LB ) 

52.26500 

DELHI (AVE) 

[BTU/LB] 

. OQOQOOOE+OO 


PRIMARY CHAMBER GEOMETRY 

ECRAT 

M 

2.682000 

RI 

[-] 

1.598000 

THETA I 

[DEGREES] 

25.42000 

RWTU 

(-) 

.4940000 

RSI 

[INCHES] 

6.900000 

RWTD 

[-1 

.2000000 

NIT 

M 

149.0000 

THE 

[DEGREES] 

8 .312000 

THETA 

[DEGREES] 

29 .22800 

EP (NOZZLE) 

[-1 

43.09743 


EXIT FLOW PROPERTIES 


P (AXIS .EXIT) 

[PSIA] 

3 . 376344 

P (WALL, EXIT) 

[PSIA] 

7 , 344903 

T (WALL. EXIT) 

(R) 

2284.013 

V (WALL. EXIT) 

[ FT/SEC j 

13825.06 

MA (WALL, EXIT) 

[-] 

4.012410 


ONE -DIMENSIONAL FLOW 

PERFORMANCE 

ISP (ODE) 

[SECONDS] 

450.0189 

ISP (ODK) 

[SECONDS] 

. 00000G0E + 00 

ISP (ODF) 

[SECONDS j 

423.3352 

CSTAR (ODE) 

[ FT/SEC] 

7392 . 1 

CSTAR (ODK) 

[FT/SEC] 

. 00000E+00 

CSTAR (ODF) 

[FT/SEC] 

7259 . 7 


TWO- DIMENSIONAL FLOW 

PERFORMANCE 

CD 

M 

.9630918 

CF (TDE) 

M 

1.914218 

CSTAR (TDE) 

(FT/SEC] 

7520.139 

THRUST (TDE) 

[ POUNDS j 

644201 . 9 

WDOT (TDE) 

[LB/SEC] 

1439.825 

ISP (TDE) 

[SECONDS] 

447.4165 


BOUNDARY LAYER PARAMETERS 


DFOPT (MABLE) 

[ POUNDS ] 

-4740, 005 

DF (MABLE) 

[ POUNDS ] 

-6740.552 

DISP (MABLE) 

[SECONDS] 

-4.548628 

THETA (EXIT) 

[ INCH] 

.3118915 

DEL* (EXIT) 

[INCH] 

.5427712 

DEL* (THROAT) 

[INCH] 

- . 6042980E 

WDOT (GG) 

[LB/SEC] 

42.06187 

EP (REGEN) 

M 

6.809522 

SQDOT (REGEN) 

[ BTU/SEC ] 

75429.73 

SQDOT (LOSS) 

[ BTU/SEC ] 

17633.74 

SUM QDOT 

[BTU/SEC] 

93063.48 

DH (SUM QDOT) 

[ BTU/LBM ] 

62.80066 


B/L PARAMETERS WITH MASS ADDITION 


DEL* BODY 
DF THETA 
DF DEL* BODY 
DF MDOT 


[INCH] 

( POUNDS ] 
[POUNDS] 
[POUNDS] 


.9672369 

13132.39 

-2000.547 

-17872.40 


THRUST CHAMBER PERFORMANCE 


THRUST (TC) (POUNDS] 

CF (TC) M 

MR (TC) (-] 

WDOT (TC) [LB/SEC] 

ISP (TC) [SECONDS] 


650942.4 

1.882642 

6.330408 

1481.887 

439.2659 
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4.0 RECOMMENDATIONS 


The focus of this work effort has been to upgrade the TDK 3 code so that it is better 
suited for use in supporting STME design and evaluation studies. Several recommendations 
are made here that would further the usefulness of the program in this respect. 

1) TDK presently requires that distributed energy release efficiency be calculated 
using input of results from another code. In practice these inputs are usually not 
available. It is recommended that a simplified method, such as that developed 
by Nurick 4 (see Appendices) be added to TDK so that a creditable procedure for 
estimating distributed energy release efficiency is available to the user. The 
method should be expanded to include all propellants and injector types of 
current and projected interest. 

2) The next logical step in developing the TDK software system is to add the 
VIPER 2 code as a module for the purpose of calculating the expansion of nozzle 
exhaust downstream of the throat plane. The need for a nozzle boundary layer 
calculation downstream of the throat plane would then be eliminated since 
viscous wall effects are included in a fully coupled manner by the PNS method. 
This would eliminate the need for a second pass through TDK. Many users of 
the code find the JANNAF procedure for core flow-boundary layer iteration to 
be confusing. Use of a PNS code eliminates this part of the procedure. SEA's 
version of VIPER currently has the capability of treating tangential slot injection 
with finite rate chemical kinetics. 

3) It is recommended that TDK be provided with a Graphics User Interface (GUI) 
for input and output. MOTIF with X- Windows should be used for this purpose 
since it is now becoming a standard for computer work stations. It appears that 
in the near future all important software application codes will be expected to 
operate in an X-Windows environment. 


T930345 


37 



REFERENCES 


1. Nickerson, G.R., Coats, D.E., Dang, A.L., Dunn, S.S., and Kehtarnavaz, H., Two- 
Dimensional Kinetics ( TDK) Nozzle Performance Computer Program , Volume I, II, and 
III, Software and Engineering Associates, Inc., NAS8-36863, dated 31 March 1989, 
prepared for George C. Marshall Space Flight Center. 

2. Kawasaki, A.H., Berker, D.R., Coats, D.E., Dunn, S.S., and Nickerson, G.R., Viscous 
Interaction Performance Evaluation Routine for Two-Phase Nozzle Flows With Finite 
Rate Chemistry, VIPER 2.0, Computer User’s Manual, Software & Engineering Asso- 
ciates, Inc., PL-TR-92-3053, January 1993 prepared for Phillips Laboratory. 

3. Nickerson, G.R., Berker, D.R., Coats, D.E., and Dunn, S.S., Two-Dimensional Kinetics 
(TDK) Nozzle Performance Computer Program Volume III, Users Manual, Software 
and Engineering Associates, Inc., NAS8-39048, dated 31 March 1993, prepared for 
George C. Marshall Space Flight Center. 

4. Nurick, W.H., DROPMIX-A PC Based Program for Rocket Engine Injector Design, 
27th JANNAF Combustion Meeting, Vol. II, pp 435-468, Cheyenne, Wyoming, 
November 1990. 

5. Nickerson, G.R., and Johnson, C.W., The Chemical Kinetics of LOXI Hydrocarbon 
Combustion, SBIR 88-1, Phase II, Contract NAS8-38468, Software & Engineering 
Associates, Inc., Carson City, Nevada, 25 July 1992. 


T930345 


38 



APPENDIX A 


Review of the DROPMIX Computer Program 


T930351 


A-l 



Introduction 


At the request of the NASA/MSFC, a review of the DROPMIX methods and computer 
program was conducted. DROPMIX is documented in Reference A-l. The computer program 
was supplied to SEA, Inc. by NASA/MSFC. It is written in the BASIC language for a 
personal computer, and was installed at SEA on a 486 PC with DOS. The program was pro- 
vided to SEA, Inc. for evaluation purposes only. A trip was taken to Irvine to discuss the 
DROPMIX analysis and its supporting data base with the author, Mr. William Nurick, who 
was most helpful and generous with his time. 

Discussion 


As discussed in Reference A-l, DROPMIX is one of a series of engineering analysis 
codes developed by W.J. Schafer Associates, Inc. for the Air Force for the purpose of assisting 
that organization in developing a booster engine for the Advanced Launch Vehicle. 
DROPMIX is the injector engineering analysis code in that series. It provides inputs to a 
thrust chamber vaporization model program, PC DER, which is based on an adaptation of the 
Standardized Distributed Energy Release program developed by Rocketdyne in 1978. 

The methods used by DROPMIX are discussed in some detail in Reference A-l, a copy 
of which is provided with this report as Appendix B. A review of the paper by SEA, Inc. 
revealed errors in several of the equations. These were shown to the author who was found to 
be aware of most of them. The copy of the paper given in Appendix B has been annotated by 
SEA, Inc. with corrections, and also with notes that are intended to clarify certain items. It 
must be stated that the errors found are nearly all typographical in that a cross-check of the 
DROPMIX source code detected no errors of any importance. 

The methods presented in the paper are empirical, and take advantage of the fact that 
measurements of spray distributions obtained using both single elements and multi-element 
injectors exhibit the same distribution whenever the measured mixing efficiency, E in , is the 
same. Thus, if E m can be predicted with reasonable accuracy, then a powerful tool is available 
for estimating combustion chamber distributed energy release efficiency. Nurick believes that 
his methods allow E m to be estimated to within -1%. However, these methods involve using 
the PC DER code for which documentation is not yet available. 


A-l. Nurick, W.H., DROPMIX-A PC Based Program for Rocket Engine Injector Design , 
27th JANNAF Combustion Meeting, Vol. II, pp 435-468, Cheyenne, Wyoming, 
November 1990 i 
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Measured mixing efficiency, E m , was originally defined by Rupe A \ and in the limit of 
an arbitrarily large number of samples corresponds to equation (1) in Reference A-l. The 
distribution function used by Nurick is given by equation (3) in Reference A-l. These func- 
tional relationships are explored in some detail in Appendix C, which also contains charts 
illustrating their behavior. 

The difficulty with an empirical approach of this type, of course, is that it cannot be 
extended beyond its existing data base with any confidence. To do so could lead to large 
errors. Furthermore, the data base is not of uniform quality. Large engines are emphasized 
over small engines, and for some injector types little or no data is presently included. 

Nevertheless, the data base contained in DROPMIX does include the results of many 
years of engine testing and the code itself is thought to provide a reasonably good description 
of the over-all mixing and atomization processes that occur in real rocket engines. The 
approach offers the possibility of becoming a highly valuable tool in the near term. That is, it 
should be applicable to the development of the next generation of large rocket engines. 

In contrast to Nurick's methods, this author is of the opinion that the CFD methods now 
being developed are not only years away from adequately modeling these processes, but also 
are very expensive to apply. The Nurick approach offers a practical way to bridge the gap that 
presently exists in this area, and will continue to exist until proven CFD analysis becomes 
available at some unknown time in the future. 


A-2. Rupe, J.H., "The Liquid Phase Mixing of a Pair of Impinging Streams", Progress 
Report No. 20-195, JPL, Pasadena, California, August 1953. 
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APPENDIX B 


27 th JANNAK Coiv.bus L i on Meeting 

Vol. II pp435-468 

Cheyenne, \JY November 1990 


DROPMIX-A PC Based Program For Rocket Engine Injector Design 


William H. Nurick 
W.J. Schafer Associates, Inc. 
Irvine, CA 


INTRODUCTION 

In 1987 the Astronautics Laboratory of Air Force Systems Command initiated development of a Personal 
Computer based engineering analysis codes for application to the Advanced Launch Vehicle engine development 
study. The codes included: (1) Engine Balance; (2) Turbopump; (3) Injector; (4) Stability; (5) Thrust Chamber 
Performance; (6) Regenerative Heat Transfer including Fins; and (7) Nozzle Performance design and analysis 
programs. The outputs from these programs have been compared against existing system component operation/ 
performance as well as proposed component designs for Step Transportation Main (STEP) Engines. These models 
have proven helpful in that independent analysis of Space Transportation Main (STEP) Engine components have 
identified several potential problem areas early in the design effort. This has allowed for design 
modifications to be made by the Contractor in a timely manner. In addition, specific deficiencies in process 
technology impacting predictive capability have also been identified. This paper describes the injector 
engineering analysis code (DROPMIX). 

DROPMIX has been developed to define propel tant spray dropsize and mixture ratio distributions produced 
by rocket engine injectors. This program utilizes the current state-of-art injector single element mixing and 
atomization correlations and includes two new relationships, one for defining spray mass/mixture ratio 
distribution and the other for defining the impact on overall injector mixing quality of intra-element mixing. 
Most injector element designs and propellant combinations of current interest are included in the program. The 
output from the DROPMIX program provides the definition of the spray characteristics for a specified injector 
design, which is input to the thrust chanter vaporization model program, PCDER. The PCDER program is an 
adaptation of the Standardized Distributed Energy Release program developed by Rocketdyne in 1978. This program 
is used to calculate overall C* performance for a given thrust chamber configuration with specified spray 
characteristics. (The PCDER program is not discussed in this paper.) The predictions from these models have 
been compared with a range of operating flight engines as well as Research/Development engines. The results 
have correlated typically within 1% with measured engine performance. The formulation of DROPMIX is described 
in this paper. 


INJECTOR SPRAY CHARACTERISTICS 

The propellant distribution and liquid dropsize(s) within the combustion chamber impacts engine 
performance, chamber material life, and combustion stability. Therefore, achievement of the required spray 
character!* st i cs is key to meeting engine mission duty cycle specifications. The propellant distribution 
characteristics are defined by injector element selection, manifold/element design, and face pattern. In this 
section the mixing and atomization processes are described relative to specific element design parameters. 
However, it is inportant to note that although the manifold design specifications can impact element mixing and 
atomization characteristics, the current model formulation does not include manifold impacts on propellant 
distributions. Consequently, the results from the DROPMIX program assume that acceptable manifold design 
criteria have been applied to the design. 

MIXING 

The interaction of the oxidizer and fuel propellant streams as they exit the injector face plane define 
the element mixing characteristics. Any geometric design or flow parameter that impacts the condition (i.e. 
shape of the exit stream, velocity, etc.) or direction of the injected propellant stream therefore impacts 
mixing. The interaction of the resulting element- to- element spray then defines the overall mixing within the 
thrust chamber. The definition of spray mixing quality is discussed below, followed by presentation of both 
single and mu l t i -element mixing design criteria and characteristics for several element types. 

Spray Mixing Quality 

Mixing quality is defined as the sum of the mass-weighted deviations in mixture ratio from the injected 
mixture ratio. This is graphically shown below (next page): 


Distribution authorized to U.S. Government agencies and their contractors; Critical Technology; November 1990. Other 
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Where: R= MR/O+MR), Injected 

r= MR/(1+MR), At any point within the spray KAfts 'fvCl-'r’) 

Note that the shaded area between the injected MR and the actual distribution curve represents the non- 
uniformity of the mixture. It is easily shown that the expression describing this deviation is: 



(i) 


Since in practice the mixing distributions are obtained by capture of finite quantities of the propellant 
simulants occupying discrete flow areas. Equation 1 can be converted to a finite difference form identical to 
that of Rupe, Ref 1 : 

MF f ( R-r, ) 

R 

r<R [ z. t r>R 

Using this expression, complete mixing is represented by a value of Em of 1 and complete non- uni f ormi ty by 0. 

Measurements of spray di str i but i ons obtained using single element and mu l t i -element injectors have been 
plotted for a wide range in overall mixing quality (i.e. Em = 50 to 90%). Typical plots are shown in Figure 
1. White all have a characteristic "s" shape the steepness of the curve, as expected, depends on the mixing 
quality, Em. Further analysis has shown that for a given Em the distribution is the same regardless of whether 
it was generated using a single or mult i -element injector; and, in addition, it is also independent of the 
element type. The spray distribution data suggested that the relatively simple expression: 

CMF - el~ b *1'!*)') ( 3 ) 



MF t ( R-r t ) 
(*- 1 ) 


( 2 ) 


satisfactorily describes the distribution. The boundary conditions (i.e. necessary to define b & f[r/R]> were 
chosen consistent with the measured distribution data which clearly shows that the curve follows a well defined 
asymptote as r approaches 1 and tends to to a? values of r %0. The resulting value for the 
constant in Equation 3 is: /S. 

YMLaA 


The functionality expression is: 


b = 0.693 

b - - XvvYu 


f (r/R) = [(r-1)/(R-1)] 

The value of n depends on the mixing quality (Em). Therefore, affixing a value to n in Equation (3) requires 
first defining the parameter Em, which relates the shape of the curve to the overall distribution. Typical 
curves for several values of Em compared with prediction using Equation 3 are shown in Figure 2. As indicated 
in Figure 2, the prediction for mixing distributions compared against cold flow data is excellent. It should 
be noted that comparisons with a wide range of element types and element configurations have been made with the 
same degree of fidelity as that shown in Figure 2. 
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Single Element Mixing 


Within the combustion environment measurement systems do not currently exist that can determine the mass 
and mixture ratio distribution produced near the injector by an element. This limitation has led to the 
development of cold flow techniques using nonreactive propellant simulants. The initial development studies 
we conducted at JPL (Ref. 2, 3, and 4), and later by the engine development Contractors (Ref 5, 6, and 7). Over 
the years the simulation techniques have been refined as it was determined that other variables such as orifice 
cavitation and gas density were important to simulate, as they impacted the spray distributions. Currently, 
measurements are being made in pressurized chambers simulating both the cavitation and chamber gas density using 
spray collection measurement systems. Unfortunately, to date there are only limited mixing data available that 
were obtained under these simulated back pressure conditions. These data, although limited, have been utilized 
to screen much of the early cold flow mixing data and reject that which is clearly invalid. The equations 
presented in this paper are considered by the Author to be the best currently available. However, it is obvious 
that these relationships should be verified as new data obtained under more representative simulation conditions 
become available and, where appropriate, changed. Currently the Astronautics Laboratory is undertaking an 
experimental program to develop mixing and atomization data for several element types. It is planned to 
incorporate the results into DROPMIX. 

Single element mixing data were used to establish several parameters: (1) the maximum attainable value of 
Em; (2) the dynamic and geometric conditions providing the maximum Em; and (3) an expression for defining the 
value of Em given the element design and operating conditions. The effort concentrated on element types and 
propellant conbinations of current interest. The mixing quality characteristics for the element types and 
propellants contained in DROPMIX are given below: 

TABLE 1 - MAXIMUM SINGLE ELEMENT MIXING QUALITY 


Propellant State 

Element Type 

Mixing Dual 

Gas/Gas 

Coaxial 

99X 

Gas/ Li quid 

Coaxial 

95X 

Like Doublet 

65X 


UnLike Triplet 

81. 5X 


Unlike Pentad 

81. 5X 

Liquid/Liquid 

Like Doublet 

SOX 

Unlike Doiblet 

sox 


Unlike Triplet 

93. 2X 


Unlike Pentad 

94. 7X 


It should be noted that these values of maximum Em represent those from specific single element configurations 
typically having orifice L/D's of 10 or greater and in some cases rounded inlets. They should therefore be 
considered only as generic values. When mixing quality data are available for a specific design it should be 
substituted for the generic value. The conditions providing maximum Em for the various elements are discussed 
below: 

Gas/Gas Propel lants 

Design criteria for mixing of gas/gas elements are limited to the coaxial element. 

Coaxial Element - A correlation for the gas/gas coaxial element describing the mixing characteristics as 
a function of geometric and flow parameters was developed by Dickerson (Ref 8) in the late 1970's. This 
correlation was successfully applied in the 1980's to advanced combustor systems using gaseous propellants. 
The correlation, however, while predicting the mixing characteristics of two concentric gaseous jets, does not 
provide a single expression for maximum mixing. Obtaining the required mixing level, therefore, requires 

a series of iterative calculations. This correlation is not currently included in the DROPMIX program, although 
there are plans to include it in the future. 

Gas/Liquid Propellants 

Design criteria for gas/liquid elements are limited to the coaxial, like impinging doublet, and the unlike 
impinging triplet and pentad elements. 

Coaxial Element - Gas/liquid coaxial elements have been divided into two types: (1) shear and (2) swirl. 
For the shear coaxial element both the gas and liquid flow axially with no tangential flow component. In the 
swirl coaxial element one or both of the propellants are swirled. Sufficient cold flow mixing characteristics 
for the swirl coaxial element are not currently available, although some limited data of Aerojet suggest that 
the mixing characterist ics for their hollow cone swirl coaxial element are insensitive to flow dynamics. 
Currently, Aerojet assumes that the mixing levels are constant at Em=80X. Other data of Rocketdyne using a 
different liquid swirler design that distributes the liquid throughout the cone, suggests that the mixing is 
similar to that of the shear coaxial element. Until more detailed correlations are developed and made available 
DROPMIX assumes that both shear and swirl elements produce the same mixing levels, regardless of the extent of 
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swirl. It should also be noted that cold flow data show that the size of the element does not significantly 
impact the element mixing quality levels. 

For the shear coaxial element the mixing characteristics from Ref 9 are shown in Figure 3. Note that the 
mixing quality was found to be a function of the gas mass flux, the ratio of liquid to gas mass flow rate, and 
the liquid velocity. The mixing level is also shown to depend on the liquid post recess. Based on these data 
the maximum mixing quality occurs at 

(^m )mj 


DR OPM IX contains both the above equation as well as a curve fit relationship for the recess of zero curve, shown 
in Figure 3. In addition, an expression relating recess to mixing at a given operating condition is also 
included in DROPMIX. 

Like Doublet Element - Like doublet gas/liquid mixing design criteria developed by Falk (Ref 10) were used 
in DROPMIX to define the impact of design parameters on mixing quality. Falk varied the fan spacing and fan 
cant angle, Figure 4. The results provided the design criteria resulting in maximcm mixing quality as: 

Fan Spacing = 0 

Fan Cant Angle =15 degrees 

DROPMIX contains relationships describing the mixing characteristics shown in Figure 4 so that single element 
mixing quality can be determined for any design condition. 

Unlike Triplet and Pentad Elements - Limited mixing data for these elements exists so the correlation 
should be considered as only interim at this time. Maximum mixing occurs when the liquid jet penetrates to the 
center of the gas jet, or: 


t 


MR V { 


^ • 

- 2000; For Recess- 1 


( 4 ) 


W o d 

(£„ @ 1.25 SIN (a) (—2*) (^) (-^) - 0.5 (5) 

Pou, 

Only limited data are available showing the inpact of off-design conditions on mixing quality for the triplet 
and none exists for the pentad. Consequently, the ability to define mixing quality for these elements in 
DROPMIX is limited. DROPMIX does, however, include off-design conditions for the triplet element. The results 
using these correlations should be used with some caution until more mixing data are available. 

L iauid/L iquid Propellants 

Design criteria for liquid/liquid elements are well doan>ented and DROPMIX includes analysis for the like 
impinging doublet element as well as the unlike impinging doublet, triplet, and pentad elements. 

Like Doublet Element - Single element mixing for the like doublet is defined as that occurring between an 
adjacent fuel and oxidizer fan, each formed by impingement of two jets of the same propellant. Extensive mixing 
studies have evaluated the impact of relative fan location, fan cant angle, and whether the fans mix on edges 
or flat sides of the fan. Studies, Ref 7 and 11, have shown that edge fan mixing provides the highest level 
of mixing quality. The correlations provided below are for edge-on-edge fan mixing only, although flat side 
fan mixing can be included at a latter date. 

Curves showing the impact of fan cant angle and fan spacing are shown in Figure 5. Note that the highest 
mixing quality occurs when the fans are in-line (i.e. 0 spacing) with an included cant angle of 45 degrees. 
Further, mixing studies have also determined the minimum flow momentum required to achieve maximum mixing. Ref 
12. These criteria are currently included in determining the maximum attainable mixing quality for an optimum 
design like doublet element. Lastly, a curve fit of the mixing characteristics shown in Figure 5 is included 
in DROPMIX for of f -opt imun design mixing level determination. 

Unlike Impinging Doublet, Triplet, end Pentad Elements - Numerous studies have been conducted to define 
mixing characteristics for unlike impinging elements. In all cases mixing quality tends to maximize according 
to a balance between the relative fuel/oxidizer momenta. Typical curves showing mixing quality characteristics 
are shown in Figure 6. The general equation that describes mixing characteristics is: 

0 ( (Ptr, I Pou,) (A^N^IA^N^) 0 E (6) 
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00.250" Fan Spacing, M = 0.34 
#0.125" Fan Spacing, M = 0.34 
A0.000" Fan Spacing, M = 0.34 



Fan Inclination Ang 




O 0.35” Spacing 
• 0.15" Spacing 
A 0.00" Spacing 



, Like Impinging Doublet Mixing Uniformity Characteristics for Several 
edge-on — edge Spacing vs. Cant Angle. 




m 


Mixing Uniformity Characteristics for Several Unlike 
Impinging Elements. 



Early work by Rupe (Ref 1), Elverun (Ref 4), and others suggested that optimum mixing occurred when: 

C=2, D=1 and E=1; Unlike Doublet, Ref 1 
C=2.75, D=1 and E=0.66, Unlike Triplet, Ref 4 
C=2.25, D=1 and E=2.75; Pentad Element, Ref 4 

The work of Elverum Ref 4 did not generate a broad data base to justify their relationships, and they noted the 
limited data upon which their correlations were based. Since that time considerably more mixing data have been 
generated (Ref 13) that has revealed that for all of the above element types mixing maximizes when: 


C=2, D=1 and E=1 


regardless of element type. This result is physically more realistic in that the same physical process governs 
mixing regardless of the number of impinging orifices. It should be further noted that regardless of orifice 
ratios, optimum mixing always occurred at E=1, although the maximum level of mixing at this condition is 
somewhat dependent on the hydraulic diameter ratio. In fact, testing over wide ranges in diameter ratio 
suggest, for circular orifices, this effect is generally small. An exception of this is for non-circular 
orifices where the hydraulic diameter ratio has a significant impact on mixing quality. In the current 
formulation, the independent impact of hydraulic diameter ratio on optimum mixing is only considered for non- 
circular orifices. Should the result of work planned by the Astronautics Laboratory (In-house) provide more 
extensive data on this variable the results will be integrated into DROPMIX. 

Equation 6 is utilized in DROPMIX to define optimum mixing for the unlike impinging elements. The mixing 
characteristic curves for the various elements are included in DROPMIX for of f- opt i rain design mixing quality 
prediction. 

Multi -Element Mixing 

The overall mixing efficiency for mult i -element burner systems is not only dependent on the single element 
mixing level but on the inter-element mixing. The degree of inter-element mixing depends on the element density 
(i.e. the spacing between the elements). A plot of overall mixing quality for several injector element designs 
ranging from single to 120 elements, as a function of the element area. density is shown in Figure 7. The 
element density is plotted in terms of injector face area per element ( inVelement ) . Note that in addition to 
optimum element designs, non-optimum elements are also shown. The straight line relationships suggest an 
Arrhenius type of equation. The general expression describing the type of relationship shown in Figure 7 is; 

E m -A In (x) ♦ ln( B) (7) 


2 

where: x = injector face area/element; in /element 

The values of A and B vary as a function of the single element mixing level (assumed to be that corresponding 
to x=100) with the resulting relationships: 




- 6.67 A ♦ 100 

(8) 


En*> 

- 2.60 B - 160 

(9) 


resulting in the final expression for the overall mixing quality of: 

. OSfe 

E m - - 1*99) ln( x ) ♦ 0.385 E^ ♦ 61.5 (10) 

Equation 10 can be used to define overall mixing uniformity given the single element mixing level and injector 
face pattern. 

It should be noted that the above correlation may not be totally descriptive and a better physical 
relationship currently being investigated is: 

E m - A ln[(l - Ato etement$(W„lA 0 )] + \n{B) (11) 

Boundary Layer Cooling 

Boundary layer cooling can be achieved in two ways: (1) the use of fuel injected near and parallel to the 
wall (i.e. film cooling) and (2) biasing the outer row of elements so that the mixture ratio differs from that 
of the other elements (i.e. core). For film cooling DROPMIX first calculates the injector distribution on the 
basis of removing the boundary layer flow from the core. Then Equation 12, adapted from Ref 14, is used to 
define the interaction of the coolant with the outer row of elements, which defines the outer zone mixture ratio 
and mass fraction. Currently the distribution of the mass/mixture ratio in the outer zone is assumed uniform, 
however the program is being modified to distribute the coolant evenly across the interaction portion of the 
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Fig. 7, Overall Injection Mixing Uniformity Characteristics 


Mixing Uniformity (E m ), % 
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outer element mass which is distributed as defined for the core. 

Mixing Limited C* Performance 

The mixing limited C* performance is defined by Equation 13. 

T| mix “ 0/C/) £ C{UR r ) dMF (13) 


Equation 13, in finite difference form, is included in DROPMIX for calculation of mixing limited performance. 
ATOMIZATION 

Atomization is known to occur in two processes: (1) primary atomization produced from interaction of the 
fuel and oxidizer propellant streams near the injector face; and (2) secondary atomization which occurs 
downstream due to drag forces applied to the primary droplet surface as the thrust chamber gas velocity relative 
to the droplet velocity increases. Both primary and secondary dropsize prediction correlations are discussed 
below. 


Primary Atomization 

The equations currently in DROPMIX defining primary dropsize for the different elements are for the most part 
based on wax studies conducted in the 60*s and 70's. The validity of the wax technique has recently been 
supported by injector element spray measurements conducted at Aerojet Techsystems and UTRC, with various 
simulants using noninterference laser diagnostic systems , Figures 8 and 9. Note the excellent correlation in 
not only the quantitative value of dropsize for the different fluids, but also the slopes; i.e. velocity and 
orifice size functionality. 

Primary dropsize equations for the various elements are given below: 


Gas/liquid 

Coaxial Element - The coaxial element is divided into two distinct types; the shear and swirl element. 
For the purposes of this report the shear element is defined by having both the fuel and oxidizer injected 
axially (i.e. without any tangential velocity component), and the swirl element is limited to only the liquid 
oxidizer having swirl. 


Shear Coaxial . The dropsize correlation for the shear gas/liquid coaxial element was taken from the work 
of Falk and Burick (Ref 15) obtained using wax/nitrogen simulants at atmospheric pressure. The data suggest 
that for low values of the correlating parameter «V -V^AMR V^] the dropsize was constant as indicated in 
Equation 14, while for larger values the dropsize decreased according to Equation 15. 




"7 


*U> /v 1 


and: 








$ 

D m - 7112 (D/ FP (14) 

whan (V - V, )I(MR V, ) < 0.55 




for. {V g -V,)!(MR V,) z 0.55 

D m - 5842 ^ \{V 0 - V, )I{MR V, ){^FP^ (1 5) 

FP - ( P/ /3) als7 (p,/.067)-°’ 4 (a,/17) 033 ( P/ /47.7)^« (16) 


* .N < 
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, Comparison of Dropsize Measurement made utilizing Malvern 
Instrument and the Correlation obtained using Wax. (p&W Ref 22) 




Swirl Coaxial. Two types of liquid swirl designs are currently employed. The first uses tangential slots 
to swirl the fluid near the entrance of the central tube producing a hollow cone atomized spray. The other 
design utilizes a mechanical swirler (i.e. helix wound metal insert) to swirl the fluid producing a full 
flowing spray cone. For the hollow cone spray the best correlation currently available is that of Aerojet (Ref 
16). The Aerojet correlation is not currently included in DROPMIX, but will be incorporated in the near future. 
DROPMIX currently utilizes the air blast atomization correlation by Lorebzetto and Lefebvre (Ref 17): 

D m - 1 .41 [(a , W, ) 0JS / ( V, p°Y p 0 , 3 J [1* W t IW g ] 17 

♦0.13 n ,[&/(*, P , ] “ [1 ♦ W, IW, ] 17 (17) 

/ 

It should be noted that there are no experimental data using swirl injection that has been used to check the 
accuracy of using Equation 17 with hollow cone swirl injector elements. Limited data do, however, suggest that 
it is a reasonable correlation for the full flow swirl injector element. 


Like Doublet Element - The like doublet element correlation included in DROPMIX is that developed by Zajac 
(Ref 18) and Nurick (Ref 19) for liquid/liquid propellants in combination with that of Zajac (Ref 20) which 
accounts for the gas velocity effects at the point of injection. The correlation is: 




> 


and for Re f 10,000 


4 

for Re 10, 


000 


ft. - t a 


* 


(18) 

. a 




D ec - 8700 V,' M D®* 7 Pj* 32 FP 




D 00 - 1.03x10* V' / ' 1 (^'& P*' PP 09) 


( 20 ) 



and final ly: 


, r 


- [1 +7xi(k)jfc v ; IL ftv, Tvgj] (21) 


Co l-2 for all cases 


, l ■ 


It should be noted that application of the above equations to define dropsize from a like doublet injector 
assumes that the gas is applied uniformly over the jet rather than an impinging fan as in current like doublet 
gas/liquid injectors. This equation should therefore be used with some caution. 


Triplet Element - Insufficient dropsize data are available for gas/liquid triplet elements. However, the 
process is not dissimilar to that of injecting a single liquid jet into a gas stream, as long as the gas jet 
diameter is infinite relative to the jet penetration length. For triplet elements where the penetration of the 
liquid jet is equal to or less than the gas jet radius the correlation of Nukiyama and Tanasawa (Ref 29) is 
utilized in DROPMIX. The Nukiyama-Tanasawi^equation is: 

n 


D m - [SQSI(V g -V f )) A + 
where 


597 B (1000 pkl(p g \))' 

A - (0,/p,) 0 - 5 
B - [p/ /(f/P / ) 03 1 0,45 


( 22 ) 


For other conditions it is assumed that the element operates similar to a like impinging element with a 
surrounding gas stream so that Equation 18 applies. The accuracy of application of these equations to the 
triplet element is highly questionable and should therefore be considered interim equations at this time. 

Pentad Element - Insufficient dropsize data exists for the gas/liquid pentad element. Since the processes 
are similar to that of the triplet element, the pentad element is treated identically to that of the triplet 
gas/liquid element. The accuracy of this approach is the same as that of the triplet and therefore should again 
be used with caution. 
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L i quid/Li quid 


tike Doublet - The like doublet element has been studied extensively by Zajac (Ref 19), Nurick (Ref 19), 
VanKleeck (Ref 21) and Hauptman (Ref 22). The resulting best fit equation, considering all of these studies, 
is that presented above as Equations 19 and 20. /y-i 

Unlike Doublet - The unlike doublet element atomization characteristics have been documented by Dickerson 
(Ref 23) and Nurick et al (Ref 13). These data are not as extensive as for the like impinging doublet and 
therefore some judgment was required in defining both the laminar/turbulent expressions as well as the jet 
development impact. The resulting equations which best fit all of these data are: 


Fuel Side 
for Re < 10000 

for Re >= 10000 


D m - 29000 V, 


',^9 P, aM D, 0 - 2 * 3 p£ ,8S {d, Id, ) 0 023 FP (23) 


D m - 34400 l',-*' 8 VV aK d,™ d 0 aoa (p,/ Po ) 04 ’ FP (24) 


Oxidizer Side 
for Re < 10000 


D m - 27200 l/^ 57 P/ s d g Pj* (d„ Id, ) * ' 7 FP (25) 


for Re >= 10000 


D m - 36500 l/,- 041 V,-° M d°* d,~° 2 ( Po /p,) a “ FP (26) 


It should be noted that these equations require further verification before they can be used with complete 
certainty. It is recommended that the user include other unlike doublet equations for comparison of prediction 
with the above expressions. Since the laser diagnostics systems can not currently discriminate between two 
fluids, it is possible that additional wax measurements can be made to better define unlike doublet atomization 
characteristics. 

Unlike Triplet - The triplet element is the least studied of the liquid/liquid elements. The atomization 
equations available are very limited in the range of parameters tested and should, at this time, be used with 
extreme caution. The equations included in DR0PM1X are: 

Fuel Side 


D m - 10400 d, a * d 0 ^ M V',- 8578 FP (27) 

Oxidizer Side 

D m - 10400 d o 0U d,- 03 * VS™ FP (28) 


The user should utilize any other equation available or hot fire data to verify the predictions before general 
use of the above equations. 

Unlike Pentad - The pentad element has been studied more extensively than the other unlike impinging 
elements although additional data would be welcomed for verification of prediction. The equations utilized for 
this element are: 

Fuel Side 


D m - 160000 d,° M d 0 -* M IV 057 V,-° M FP (29) 
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Oxidizer Side 


* 

D m - 36000 d, 0 ' d° Ai ^ (oH) V," 089 FP (30) 


Secondary Atomization 

The secondary atomization correlations are also based on wax measurements; however, these correlations are 
lacking in several key physical process descriptions. For example, the correlations are only in terms of 
input/output parameters correlated against the maximum gas velocity achieved in the test chamber. Application 
of the secondary breakup equation requires selection of a "maximum 11 gas velocity for breakup and, further, that 
all of the breakup occurs in less that 2 inches from the injection plane. The latter assumption appears to be 
experimentally verified. For small space engines having chamber lengths of 1 to 4 inches, the exact location/or 
conditions producing secondary breakup are required. The selection of the maximum gas velocity is also critical 
to the degree of secondary breakup achieved. However, insufficient data are available to discern whether this 
velocity should be that consistent with the maximum delta-V between the droplet and gas, that occurring in the 
charter at the start of convergence, or that at the nozzle throat. Another complication is that in cold flow 
tests vaporization of the droplet does not occur. Consequently, the only parameter impacting breakup is the 
dynamic interaction. However, in the combustion chamber the dropsize is changing due to both dynamics and 
vaporization. 

Inspection of the secondary breakup characteristics shown in Ref. 25, suggests that the secondary dropsize 
decreases rapidly as the gas velocity is increased and then reaches a minimun dropsize. The gas velocity 
corresponding to this minimum dropsize is reasonably small (i.e. about 300 ft/sec) suggesting that for normal 
large engines with contraction ratios in the range of 2 to 3, secondary breakup will occur quickly and reach 
the minimum secondary dropsize consistent with the initial primary dropsize. DROPMIX currently utilizes the 
gas velocity at the start of convergence (i.e assuming complete vaporization) which assures that the minimum 
secondary dropsize is always attained. An option is also included that allows for "short" chamber lengths 
(i.e. less than about 6 inches), for secondary atomization to occur in two steps, first in the chamber and 
secondly in the nozzle where the maximum gas velocity at the throat is used for the second breakup calculation. 
To be physically realistic in cases where breakup can occur in two steps, PCDER should be run to the start of 
convergence then the breakup expression applied to the remaining spray before proceeding into the convergence 
section of the nozzle. For small space engines where the total length to the throat is about 2 inches or less 
secondary breakup is applied in one step using the nozzle throat velocity as the maximum gas velocity in the 
secondary breakup expressions. 

Two different breakup equations are included in DROPMIX. The first correlation is for the coaxial element 
and the other for all other elements. The shear coaxial element secondary droplet breakup characteristics, 
developed by Falk (Ref 24) are shown in Figure 10. The resulting best fit equation of the correlation line is 
given below: 


D, ID p - 0.2314+0.3185 S-0.0461 S 2 +0.0016 S 3 (31) 

where S - ln[(0 p /110) 3 ABS ( V g - V , )] (32) 


The secondary breakup equation for the other elements is that of Zajac (Ref 25): 

for -1 s A VIV f s 1.25 

/, - £^[ 1 - 1 . 77 x 70- 3 £|p LVIV, e ( " M44V/V ' ,UV/V ''') ( 33 ) 

? T V vr 

for &V IV, > 1.25 J- 

* a 

D, - t £^[1-1.52*70- 3 Dfe ] - 12 ln(AV/V,) (34) 

1 I 

? T 

MANIFOLD/ORIFICE DISCHARGE COEFFICIENT 

The manifold/orifice discharge coefficient can vary significantly from one design to another. For most 
impinging element injectors the propellant enters the manifold either: (1) directly at the rear of the 
injector, is diffused into a dome like manifold then distributed to the individual element orifices via axial 
step drilled passages, or (2) into a ring manifold surrounding the injector diameter, then fed directly into 
radial step drilled passages to the individual element orifices. In general these element orifice inlets are 
sharp edged, with relatively short L/D‘s (i.e. 2 to 5). It should be noted that the passages intersect the 
element orifices at some angle (i.e. 30 to 120 degrees are typical). For these types of designs the discharge 
coefficient at typical engine chamber pressures will vary from about 0.5 to 0.85 depending on the specific 
design conditions. Coaxial element designs are more complicated with control orifices or tangential slots at 
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the liquid entrance to the central post as well as slots on the annulus of the gaseous outer sleeve. For these 
designs the overall discharge coefficient can range from 0.2 to 0.9. The myriad of possible designs makes 
incorporating a generalized manifold C d routine impractical. 

DROPMIX includes equations for calculating discharge coefficients for sharp edged orifices as a function 
of L/D and Reynolds Nunber, in addition to the liquid post and gaseous annulus discharge coefficient for coaxial 
elements. The manifold pressure drop is input by the user. The user also has the option to input any value 
of orifice discharge coefficient for his particular design. 


The orifice discharge coefficient is defined from Figure 11, after Lichtarawitz (Ref 25) with the resulting 
equation: 

MC d - 1 1C d + 20 (Re h (1 +2.25 LJd) - 0.005 L 0 Id I A (35) 


«v- '%*■ 

- -/» 

zOJ-- \j- ^ 

-, v - 

- l 


where A - [1+7.5 ( k>ojjo0015 Re h ) 2 ] 

Re, - [2 g p A/’] os d (36) 

for 2 s L„ Id s 10 ; 10 s Re, s 20,000 


- 0.827-0.0085 L„ Id 


(37) 






For Re > 2x104, C d is equal to C d(J . This equation is limited to L 0 /d of 2 or greater. Calculations of C d from 
Equation 31 have been verified against data taken with typical injector designs. 

The procedure for determination of discharge coefficient for coaxial elements is also included in DROPMIX. 

ANALYSIS PROCEDURE 


The above equations are programed into an analysis procedure which can be used to define the "optimum" 

design, off-optimum design, or simply predict the mixing quality and dropsize for an existing or proposed 

injector design. An overall flow diagram for DROPMIX is provided in Figure 12. 

The initial calculation defines the ''optimum" injector design for the given input requirements (i.e. 
element type, propellant, engine operating conditions, and overall mixing quality). In this initial calculation 
the minimun delta-p is set at 20% of chamber pressure based on chug stability criteria. The other propellant 
delta-p is defined from optimtm mixing requirements. The number of elements per square inch of injector face 
area is determined from the mixing quality limit selected. The injector can include film cooling or outer zone 
mixture ratio bias. Currently the interaction of the outer boundary flow with the core flow, to determine the 
overall mixing quality as well as the mean gas temperature near the wall, is defined by the % interaction. Once 
the mixing quality, distribution, and element sizes are determined, the resulting primary and secondary 
dropsizes are specified for an optimum element design. 

After the initial calculat iona l procedure is completed then the user can change a number of variables 

sequentially to either investigate the impact of the changes on mixing and atomization or to represent an actual 

injector design. The variables that can be changed are: 

o C. 

o Element Type 
o Mixing Level 
o Delta-p 
o Thrust /Element 
o Chaober Gas Velocity 
o BLC Parameters 


When the user is satisfied with the specifications, an output is requested. A Typical output is given in 
Table 1. Note that in addition to the input specifications, all injector design parameters, dimensions, and 
dropsizes are also given. Of particular interest for input into PCDER are the values of secondary dropsize and 
the mixing parameter "n" defined in Equation 3. In the event that boundary layer cooling or mixture ratio bias 
is included in the design, then the various mixture ratios and mass fraction in the boundary layer are also 
required inputs to PCDER. 


MODEL FIDELITY 

A verification of the ability of the DROPMIX Engineering Analysis program to adequately predict the 
injector mixing quality and atomization was accomplished by utilizing the DROPMIX output as input into PCDER, 
calculating the energy release efficiency then comparing the predictions against actual engine data. To be 
thorough, flight as well as research engines were analyzed. The results are described below. Currently, small 
space engines have not been analyzed in sufficient detail to verify the fidelity of the formulation for this 
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class of engines. It is planned that this will be accomplished in the next year. 


FLIGHT ENGINES 

Several large launch flight vehicle main engines were analyzed to compare measured performance with that 
predicted by entering DROPMIX specifications into PCOER. The engines analyzed are: 

o J-2 (LOX/HO 
o J-2S (L0X/H 2 ) 
o SSHE (LCDC/HO 
o F-1 (LOX/RP-1) 


In each case the injector and thrust chamber design as well as operating conditions were obtained and input 
into DROPMIX to define both the mixing uniformity and liquid propellant dropsize. The output from dropmix was 
then input into PCDER and overall mixing efficiency determined. For the LOX/H 2 engines all injector designs 
were shear concentric tube elements. The LOX/RP-1 engine utilized like-impinging doublet elements. The 
analysis for all of these injectors revealed that vaporization was virtually complete within the combustion 
chamber requiring that the resulting performance is that defined by mixing. The results of the analysis are 

” _ . . . ^ . i; .i. I i ^ i i i » k i n It thaf (tlcj a c i iron T h AC A 



adequately predict mixing limited efficiency. 


RESEARCH ENGINES 

Several research injector/thrust charrber assemblies were evaluated wherein both mixing and vaporization, 
limited the C* performance. Those evaluated include: 

o MSFC (LQX/ CH,) 
o Aerojet (LOX/RP-1) 
o Rocketdyne (NTO/50-50) 
o Rocketdyne (LOX/RP-1) 
o Rocketdyne (L0X/RP-1/CH 4 ) 

MSFC L0X/CH 4 Study 

The NASA-MSFC LOX/CH, TCA design evaluated utilized several shear concentric tube injectors, described in 
detail in Ref 27. Tests were conducted over a range of about 2 to 4 in mixture ratio and 1200 to 2400 psia in 
chamber pressure. For these injector/ thrust charter designs DROPMIX/PCDER predicted 100% vaporization 
efficiency as in the flight engines discussed above. It should be noted, however that these test engines 
utilized a different fuel propellant. Consequently, as for the flight engines these predictions are based only 
on the mixing limit as impacted by element design, mixture ratio, and chamber pressure. The results are shown 
in Figure 13. Note that the comparison between measured and predicted C* are in general excellent. Also note 
that the two points that do not compare within 1/2 % are clearly inconsistent with the other measured C 
efficiencies and are therefore suspect. 


Aerojet LOX/RP-1 Study 

A detailed cold flow/hot fire test program was conducted at Aerojet, Ref 28, as part of the Advanced Launch 
System Advanced Development Program studies. In these studies both unlike impinging triplet and like impinging 
doublet injectors were tested. Tests were conducted over a wide range in chamber length ( 8 to 20 inches). 
Mixture ratio and chamber pressure variations were relatively small. In the analysis the actual cold flow 
measured mixing uniformities <i.e. Em) were input to DROPMIX. A comparison of the predicted and measured 
performance are provided in Figure 14. Note that the predictions are within IX of the measured values. Also 
note that the C* predictions with chamber length virtually match the hot fire measured values suggesting that 
the vaporization formulation in PCDER is technically acceptable. 

Rocketdyne NTO/50-50 Study 

This study (Ref. 30) conducted in 1967 tested a number of like doublet, unlike doublet and pentad elements 
firing NTO/50-50 propellants at 5000 LBf . Tests were conducted at 300 psia chamber pressure and chamber 
lengths <L') of between about 4 and 16 inches. For the unlike impinging elements, blowapart was experienced 
which makes that data unacceptable for analysis. The analysis was therefore limited to the like impinging 
elements. It should be noted that blowapart criteria has been developed and can be included in Dropmix. 

Both six element and 100 element configurations were tested. Three configurations of the six element 
injector were studied which were identical except that the element fan-to-fan spacing was varied from zero (i.e. 
in-line) to 0.35 inch fan spacing between adjacent fuel and oxidizer fans. The 100 element configuration was 
designed with zero fan spacing. 

The six element injectors were tested with chamber to throat lengths (L 1 ) of 4.57, 9.57, 13.97, and 19.97 
inches. A comparison between the predicted and measured C* performance is given in F [ 9 ure Note the 
excellent agreement for all the data except for the two lowest performance values (i.e. about 65% C ). Based 
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TABLE 1 . TYPICAL DROPMIX OUTPUT 


ENGINE SPECIFICATIONS 


CHAMBER PRESSURE, psia 1974 

SEA LEVEL THRUST, lb 2631 3 

MIXTURE RATIO 3.52 

OXIDIZER FLOWRATE, Ib/sec 68.78 

FUEL FLOWRATE, Ib/sec 1 9.54 

THROAT DIAMETER, in 3.31 

CHAMBER DIAMETER, in 5.66 

CHAMBER GAS VELOCITY, ft/sec 881 
CHAMBER GAS DENSITY, Ib/ft3 0.57 

SPECIFIC HEAT 1.24 


SPRAY SPECIFICATION 
LOX: 

VELOCITY, ft/sec 

133 

PRIM. DROPSIZE, micron 

94 

SEC. DROPSIZE, micron 

77 

FUEL: 

VELOCITY, ft/sec 

614 

PRIM. DROPSIZE, micron 

0 

SEC. DROPSIZE, micron 

0 

MIXING EFFICIENCY, % 

98.73 

DISTRIBUTION FACTOR 

6.5001 


M/MR DISTRIBUTION 

MMR MMF 
2.156 0.0058 

2.341 0.0248 

2.549 0.0656 

2.786 0.1187 

3.059 0.1598 

3.374 0.1706 

3.745 0.1519 

4.186 0.1173 

4.720 0.0806 

5.380 0.0502 

6.216 0.0285 

7.309 0.0147 

8.799 0.0068 
10.951 0.0027 
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INJECTOR SPECIFICATIONS 


ELEMENT TYPE: COAXIAL LIQUID/GAS 


NUMBER OF ELEMENTS 

82 

THRUST PER ELEMENT, lb 

321 

OVERALL MIXING, % 

90 

SINGLE ELEMENT MIXING, % 

51.9 

LOX: 

POST WALL THICKNESS, in 

0.035 

POST ID, in 

0.134 

POST OD, in 

0.204 

POST RECESS/DL 

1.00 

Cd 

0.47 

DELTA-P, Ib/in2 

560 

FUEL: 

GAP, in 

0.010 

Cd 

0.7 

DELTA-P, Ib/in2 

684 


PROPELLANT SPECIFICATIONS 


PROPELLANT COMBINATION: 

LOX/CH4 

OXIDIZER: 

TEMPERATURE. F 

-249 

DENSITY, Ib/ft3 

64.78 

VISCOSITY, cp 

0.12 

FUEL: 

TEMPERATURE, F 

24 

DENSITY, Ib/ft3 

8.24 

VISCOSITY, cp 

GAS 
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TABLE 2. COMPARISON OF DROPMIX C* PREDICTION VS. 
FULL SCALE ENGINE MEASURED C. 





























on the DROPMIX predictions, except for these two tests, all of the measured input values compare. For the two 
low performance tests a comparison between delta-p and flowrate was not obtained, suggesting that there is some 
inconsistency between the measured input data. This overall comparison is extremely encouraging since 
performance varied from 90% down to 65X resulting in an excellent check on the ability of DROPMIX to handle both 
mixing and vaporization limited combustion, as well as the validity of the mixing and atomization correlations 
(i.e. for like impinging doublets). 

The 100 element injector was tested at only two chamber to throat lengths, 9.57 and 4.57 inches. Only the 
results of the 9.57 inch testing is presented since the data at 4.57 did not provide a reasonable check on the 
validity of the input data. Note that for the 9.57 inch chamber, the predicted and measured C* performance also 
check within 1%. 

Rocketdyne LOX/RP-1 Study 

During this study (Ref. 31) a total of five different injector designs were tested. Of these only the 16 
element O-F-O triplet and 60 element circumferential fan tike doublet injectors could be evaluated using built- 
in DROPMIX mixing correlations. The other designs utilized element designs and configurations not currently 
included in the program. It should be noted that, since dropsize correlations for most of the configurations 
are included, predictions could have been made if single element mixing experiments had been accomplished. The 
nominal test conditions are 2000 psia and L* of 19 inches. 

The results for the triplet and like doublet are shown in Figure 16. Note that there are two data points 
for each injector. The open symbol are test results from earlier nominal 2 second steady state run testing and 
the solid symbol from a half second steady state. The like impinging doublet injector test results were nearly 
the same for both test series while there was a substantial difference between the two runs for the triplet 
injector. A reasonable explanation for this difference is not known at this time. Additional tests using the 
same hardware obtained in another study will be analyzed to determine consistence of the test results. 

Rocketdyne LOX/RP-1/CH 4 Study 

This study (Ref. 32) utilized two propellant combinations (LOX/RP-1 ,CH,) at chamber pressures of 2000 psia, 
chamber diameter of 3.5 inches and L' of 19 inches. Only the LOX/CH^ results will be discussed in this paper, 
however, it is planned that the LOX/RP-1 tests results will be analyzed at a later date. The LOX/CH, injector 
was a 58 element concentric tube injector. Three configurations of this basic injector were utilized: (1) 
uniform elements; (2) mixture ratio bias of the outer ring of elements and; (3) uniform injector elements with 
film cooling injected at the beginning of convergence. The mixture ratio bias configuration can not be analyzed 
at this time since DROPMIX can not currently handle two separate distributions. 

A comparison of the predicted and measured test results are shown in Figure 17. Note that both the 
baseline and film cooling results compare reasonably well with the predictions. 

CONCLUSIONS 

The results are encouraging in that the basic engineering formulation appears to be a reasonably good 
description of the mixing and atomization processes. While comparisons have been made using a wide range in 
engine designs, injector thrust/elements, and propellant combinations, not all of the element type correlations 
have been verified. As noted in the body of this paper, many of the correlations are lacking in an extensive 
data base and processes such as secondary breakup need further development. A need also exists to demonstrate 
the ability of the model to adequately predict mixing/atomization characteristics for small engine applications. 

The Air Force Astronautics Laboratory is currently planning to conduct an extensive cold flow study in 
their new cold flow facility at the Astronautics Laboratory. The facility is capable of pressurized atmospheres 
of up to 2000 psia, and is equipped to measure both mixing uniformity and atomization. An additional study of 
secondary breakup is also planned. The correlations developed from these studies will be incorporated into 
DROPMIX and should provide most of the needed data. A large part of the study will also concentrate on 
developing manifold design criteria to ensure that these impacts can be accounted for in DROPMIX. 


Lastly, the program has been loaned to several users in an effort to obtain independent checkout of the 
predictions. Feedback has been favorable, but application has been limited to large engines. 
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Fig. 14, Comparison of Aerojet Data vs. Dropmix/PCDER Prediction. 







Fig. 15, Comparison of Rocketdyne Research Engine Data vs. Dropmix/PCDER Prediction 
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Fig. 16, Comparison of Predicted and Measured c* Efficiency for two 
Different Element Types. (Rocketdyne) 





oo 

cn 


Fig. 17, Comparison of Measured and Predicted c* Efficiency for 
a 58 Element Concentric Tube Injector. (Rocketdyne) 
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NOMENCLATURE 


A = Orifice Area 

b * constant (0.693) 

* 

C = C-Star 

C^ a Discharge Coefficient 

CMF = Cumulative Mass Fraction 

d = Orifice Diameter 

D = Drops ize 

E m = Overall Mixing Uniformity Defined in EQ. 1 

E mo = Sl * n 9l e Element Mixing Uniformity 

FP = Physical Property Correction defined in EQ. 16 

L = Length 

L' = Distance from Injector Face to the Throat 

MF = Mass Fraction 

MR = Mixture Ratio (W o /W f ) 

N = Number of Fuel or Oxidizer Orifices in Element 

Pq = Dynamic Pressure Ratio 

Pj = Velocity Profile Correction Factor - Ratio of Centerline to Mean Dynamic Pressure 

A P = Orifice Pressure Drop 

r = MR/(1+MR) # Local Value 

R = MR/O+MR), Injected Value 

Recess = Oxidizer Post Setback from Fuel S 
Re = Reynolds N ember 

V = Velocity 

W = Flow Rate 

Fan Spacing - Spacing Between Adjacent Fuel and Oxidizer Edges 

Fan Cant Angle - Total Induced Angle Between Adjacent Fans relative to the Injector Face 


Greek 


a = Orifice Element Impingement Half Angle 

* 

*1^ - Mixing Limited C Efficiency 


P = Viscosity 

p = Density 

o = Surface Tension 


or 




Subscripts 
c = Chamber 

co = Primary with Zero Gas Velocity 

f = Fuel 

9 = Gas 

H = Hydraulic 

inn = Inner 

J = Jet 

l 55 Liquid 

m = Mass Mean 

max = Maximum Value 

o = Oxidizer 

out = Outer 
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p = Primary 

r = Relative 

s = Secondary 

t = Theoretical 

TC = Total Core 
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Rocket Engine Injector Spray Mixing Quality 
File 

G. R. Nickerson and C. W. Johnson 
5 February 1992 


ABSTRACT 


A method for predicting propellant spray drop size and mixture ratio distributions 
produced by rocket engine injectors is presented by Nurick in Reference 1. The method is 
empirical, and takes advantage of the fact that measurements of spray distributions obtained 
using both single element and multi-element injectors exhibit the same distribution whenever 
the measured mixing efficiency, E m , is the same. The purpose of this memo is to discuss the 
properties of the distribution function used by Nurick to characterize the sprays. The 
relationship between the distribution function and the Rupe mixing efficiency, E m , is also 
discussed. Charts of E m as a function of the oxidizer fraction for the total propellant are 
presented for a series of shape factors. 

Liquid Oxygen (LOX)/Hydrogen propulsion systems are normally operated at a mixture 
ratio of approximately one for fuel rich gas generators and pre-burners, and at a mixture ratio 
of approximately six for the main combustion chamber. For these operating conditions, charts 
are given that show both the mass density function, and the mass distribution function for a 
series of shape factors. 


Ref. 1. Nurick, W. H., "DROPMIX - A PC Based Program for Rocket Engine Injector 

Design," 27th JANNAF Combustion Meeting, Vol. n, pp 435-468, Cheyenne, 
WY, November 1990. 
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NOMENCLATURE 


CMF 

dm/dr 

E m 

f(r) 

m 

m i 

m 2 

MFi 


MR 

n 

r 

R 

r 

max 

density 


Cumulative Mass Fraction for the distribution. Jo mf dm = f(r) 

Mass Density function, dm/dr = df(r)/dr 
Rupe mixing efficiency, see eqs (1) and (5) 

Cumulative Mass Fraction, m = CMF 

Equivalent area displacement height, r < R. See Fig. 4 

Equivalent area displacement height, r > R. See Fig. 4 

Fraction of mass measured in compartment i relative to the total mass 
measured. See eq. (1) 

Engine mixture ratio. MR = R/(l - R), 0 < MR < » 

Shape factor exponent for f(r) 

oxidizer fraction, 0 < r < 1 

Engine over-all oxidizer fraction. R = MR/(1 + MR), 0 < R < 1 

r value within the distribution where dm/dr is maximum, and where f(r) has 
inflection. See eq. (4) 


Mass Distribution function. f(r) = exp 


H 


T920203 


C-2 



SPRAY MIXING QUALITY 


Cold flow data is often available for rocket engine injectors so that data on 
injector/thrust chamber mixing efficiency is available. Typically, the simulants for the 
propellant oxidizer and fuel are flowed through the injector and collected downstream in a 
compartmentalized container. The oxidizer and fuel content of each compartment is 

determined, and a mixing efficiency, E m , is calculated by the following formula: 

E m = 1 - i I MF;(R - r,) - ^ I MFj(rj - R) (D 

i= 1 j=l 

r j <R r j >R 


where 


R 

r j 

MFj 

MF. 


is the oxidizer fraction for the total flow 
is the oxidizer fraction for the i th compartment, and q < R 
is the oxidizer fraction for the j th compartment, and rj > R 
is the mass fraction found in compartment i 
is the mass fraction found in compartment j 


The total number of compartments is 


n t = n 1+ n 2 


The mixing efficiency, E m , given above is due to Rupe 2 , and is known as "the Rupe 
mixing efficiency.” 


Ref. 2. Rupe, J. H., "The Liquid Phase Mixing of a Pair of Impinging Streams," 

Progress Report No. 20-195, JPL, Pasadena, CA, August 1953. 
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The engine mixture ratio, MR, and engine oxidizer fraction, R, are related as 


R = MR/(1 + MR) 0 < R < 1 

MR = R/(l -R) 0 < MR < » 


The distribution function for the spray is an expression that gives cumulative mass 
fraction, m, as a function of oxidizer fraction, r. The distribution function, f(r), given by 
Nurick 1 contains two parameters, R and n, and is 


CMF 

J 


dm = e^ n ^ 


o 



= f(r) 


(2) 


The mass density function, (dm/dr), for this function is found as follows: 


I 

o 


dm 

dr 


dr = f(r) 


d } 

dm ' 

dr = 

dm 

0 

h d r 


^ d r 


dr = df(r) 


and 


dm ' _ df(r) 
dr _ - dr 


' dm ' 

_ -n In ^ 

r 1-r ' 

_ dr ^ 

1 - R 

1-R 


f(r) 


(3) 


Example plots for equations (3) and (2) are shown in Figures 1 and 2, respectively. 
Note that for r=R, equation (2) gives a cumulative mass fraction (CMF) value of 1/2, i.e., in 
Figure 1 the area under the curve to the left of R=r is equal to the area under the curve to the 
right of R=r. Note also that Figure 2 has a total area of one. The distribution function is 
always greater than zero at r=0, equal to 1/2 at R=r, and equal to one at r=l (all mass 
accounted for). 
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dm 

dr 


a lhs = a rhs 

-L 


too im . ioo -*ac '-oo *oo .i{« 
oxidizer fraction, r 


•00 t OOO 


Figure 1. Mass Density Function, (dm/dr) = - n In \ (1-R) n (l-r) n 1 f(r) 
(In this example R =.5, and n = 6.22485) 



Figure 2. 


CMF 

Distribution Function, m = f(r) = I dm = exp 

o 

(In this example R =.5, and n = 6.22485) 



d CMF 
dr 


CMF 

R)4 m=.5 

} 
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For the distribution function to have the characteristic "S" shape that is observed in 
spray distribution measurements, it is necessary that the mass density function have a 
maximum value on the interval 0 < r < 1. The maximum are at 



and it is found that 


dm 
. dr 


max 


is at 


r 

max 

density 


= 1 - (l-R) 


r -1 

l in - ? 



(4) 


For this value to be real, it is required that n be greater than one. Equation (4) gives 
the inflection point seen in Figure 2. The inflection point occurs at r = R only when nn»or 
at 


n = l/( 1 + In \ ) = 3.25889 


The minimum value for the inflection point is less than R, and only occurs when 


(n - 1) In 


-1 n-1 

In T n 


= 1 


which is at 


n = 6.22485 

The corresponding maximum density value for r, 

= 1.03122R - 

min 


” r 

= 


max 


f(x) 

density 

min 

inflection 


is 

.03122 
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As the mixing efficiency approaches unity, the areas Aj and A 2 shown in Figure 3, 
below, approach zero. 


m 



Figure 3. The Areas A! and A 2 that are Used to Indicate Mixing Efficiency 


The hatched area shown in Figure 3 is 


R 

A = Aj + A 2 = I f(r) dr + 
o 


(1-R) - 1 f(r) dr 

R 


The Rupe mixing efficiency, E m , is defined in terms of the displacement heights, mj 
and m 2 , as shown in Figure 4. 



Figure 4. Displacement Heights, mj and m 2 , Used to Define the Rupe Mixing 

Efficiency, i.e., E m = l-mi-m 2 
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where 


rri! R = Aj and m 2 (1-R) = A 2 . 

The geometric definition for E m is 

E m s 1 - mi - m 2 

so that 

E — 1 - 
m R 1-R 

R 

Aj = J f(r) dr 

o 

1 

A 2 = (1-R) - J f(r) dr 

R 

R r 1 1 

E m = 1 - 1 I f(r) dr - 1 - R - I f(r) dr 

0 R 

i 1 i R 

E m = y~^ } f(r) dr - i J f(r) dr (5) 

R o 

The above expression is equivalent to equation (1) in the limit when an arbitrarily large 
number of samples are taken. It is also equivalent to equation (1) of Reference 1. 

If the Rupe mixing efficiency, E m , for a spray is known, then the spray distribution can 
be obtained using equation (5). This is done by setting R = MR/(1+MR) and finding a value 
of n that satisfies equation (5). Once n is known, the spray distribution, f(r), given by equation 
(2) is known. Equation (5) is solved by numerically integrating f(r), and iterating values of n 
until the Right Hand Side equals E m . Plots of E m vs R for various values of n are presented in 
Figures 5 and 6. These charts can be used to conveniently estimate a value for n from given 
values of E m and R. 
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The mass density function, eq. (3), is plotted in Figures 7 and 8 for a mixture ratio of 
one, i.e., R = 1/2 (typical of fuel rich O 2 /H 2 gas generators). Figure 7 shows the function for 
values of n < 8, and Figure 8 shows the function for values of n > 8. The maximum curve for 
the function is also shown, as calculated by equation (4). Oxygen/Hydrogen systems are of 
primary interest, and for these systems the combustion chamber is usually operated at a 
mixture ratio of approximately six, i.e., R = .857. The mass density function for MR=6 is 
given in Figure 9 for values of n > 8. 

Figures 10, 11, and 12 show the distribution function, f(r), corresponding to Figures 7, 
8, and 9, respectively. 

As mentioned earlier, the mass distribution function, f(r), is always greater than zero at 
r=0, i.e.. 


f(0) = CMF(O) = exp { [In (1/(1-R)) n } 

In Figure 13, CMF(O) is plotted vs R for several values of n. Note that 

lim CMF(O) = 1/2. 

R- 4 O 
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Figure 5. Rupe Mixing Efficiencies for Different Shape Factor Exponents (n). 
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Figure 7. Mass Density Function vs r for n < 7, R=0.5 
(i.e. MR=1, typical of Fuel Rich LOX/H2 
Gas Generators and Preburners) 




locus of dm/dr max. 



Figure 8. Mass Density Function vs r for n > 8, R=0.5 
(i.e. MR = 1, typical fo Fuel Rich LOX/H2 
Gas Generators and Preburners) 
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Figure 9. Mass Density Function vs r for n 2 8, R-0.857 

(i.e. MR = 6, typical to LOX/H2 
Primary Combustion Chambers) 
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Figure 10. Mass Distribution Function vs r for n < 8, R=0.5 
(i.e. MR = 1, typical of Fuel Rich LOX/H2 
Gas Generators and Preburners) 
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Figure 1 1 . Mass Distribution Function vs r for n > 8, R = 0.5 
(i.e. MR = 1, typical of Fuel Rich LOX/H2 
Gas Generators and Preburners) 
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